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Abstract 

This paper addresses the issue of integrable structure in a modi- 
hed mehing crystal model of topological string theory on the resolved 
conifold. The partition function can be expressed as the vacuum ex- 
pectation value of an operator on the Fock space of 2D complex free 
fermion helds. The quantum torus algebra of fermion bilinears behind 
this expression is shown to have an extended set of "shift symmetries" . 
They are used to prove that the partition function (deformed by exter- 
nal potentials) is essentially a tau function of the 2D Toda hierarchy. 
This special solution of the 2D Toda hierarchy can be characterized by 
a factorization problem of Z x Z matrices as well. The associated Lax 
operators turn out to be quotients of hrst order difference operators. 
This implies that the solution of the 2D Toda hierarchy in question 
is actually a solution of the Ablowitz-Ladik (equivalently, relativistic 
Toda) hierarchy. As a byproduct, the shift symmetries are shown to 
be related to matrix- valued quantum dilogarithmic functions. 
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1 Introduction 



This is a sequel of our previous work [1] |2] on the integrable structure of 
the "melting crystal model" of topological string theory [3] and 5D M = 1 
supersymmetric U{1) gauge theory [3]. We here address the same issue for 
a modified melting crystal model that is related to topological string theory 
on the resolved conifold. In particular, we present an affirmative answer to 
a conjecture in our preliminary study on this issue |5]. In the perspectives 
of topological string theory, these two models are the simplest non-trivial 
cases of the topological vertex construction [6] of string amplitudes on non- 
compact Calabi-Yau 3-folds. The Calabi-Yau 3-folds relevant to our models 
are the first two members of the so called "local CP^ geometries" Xi = 
0{-l) © 0{l - 2) -> C¥\ I = 0, 1, . . .. Local Gromov-Witten invariants 
of these manifolds are studied by the localization technique Moreover, 
the genus- zero part of those invariants is determine by the random matrix 
technique [H]. 

The partition function Z of the previous model, see fl2.1l) . is a sum of 
weights labelled by integer partitions of arbitrary lengths. This statistical 
sum is derived from the partition function of random plane partitions (3D 
Young diagrams) by the method of diagonal slicing [9|. The aforementioned 
integrable structure [HIS] emerges when Z is deformed by external potentials. 
The partition function of the deformed model is a function Z{s, t) of a discrete 
variable s and a set of continuous variables t = (ti,t2, • • •)• -5 is a parameter 
of the external potentials and t^s are coupling constants. The deformed 
partition function Z{s,t) turns out to be related to a tau function r(s,t) 
of the ID reduction of the 2D Toda hierarchy [TOl [11]. s and t play the 
role of a lattice coordinate and time variables in the ID Toda hierarchy. 
The relation to the Toda hierarchy is explained with the aid of 2D complex 
free fermion fields that are extensively used for integrable hierarchies as well 
[121 US]- Z{s,t) can be thereby expressed as the vacuum expectation value 
of an operator in the fermion Fock space. We found that this fermionic 
expression of Z{s, t) can be converted to the tau function r(s, i) by a change 
of variables and a simple multiplicative factor. A technical clue therein is the 
notion of "shift symmetries" [Jj, i2j generated by the adjoint action of special 
"vertex operators" [9j on the quantum torus algebra of fermion bilinears. 

Our approach to the partition function Z' of the modified model, see 
fl2.7l) . is mostly parallel, but has a new feature as well. Namely, Z' is de- 
formed by two sets, rather than a single set, of external potentials. The 
deformed partition function is a function Z'{s,t,i) of a discrete parameter 
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s and coupling constants t = (^1,^2, • • •) and i = {ii,t2, . . .) of the exter- 
nal potentials. We introduce an extended set of shift symmetries (related to 
vertex operators of another type [H]), and show that the deformed partition 
function Z'{s, t, i) is related to a tau function r'(s, t, i) of the 2D Toda hier- 
archy in the same way as Z{s, t) is related to r(s, t). It is natural to expect 
that this tau function, too, is a solution of a reduced system of the 2D Toda 
hierarchy. This is indeed the case. As we conjectured in the preliminary 
study [5], we here prove that r'(s,t,t) is a solution of the Ablowitz-Ladik 
hierarchy [15] or, speaking more fairly, the relativistic Toda hierarchy [T6] . 
This explains an origin of the integrable structure that Brini [17j observed 
in the generating function of local Gromov-Witten invariants of the resolved 
conifold by a genus-by-genus analysis. 

The Ablowitz-Ladik hierarchy is a spatial discretization of the nonlin- 
ear Schrodinger hierarchy. A traditional Lax formalism is based on zero- 
curvature equations of 2 x 2 matrices \15\. A bilinear formalism consists 
of Hirota equations for three tau functions r, a and a [18] . In the context 
of random matrices, a semi-infinite lattice version emerges in the unitary 
matrix model and the associated biorthogonal Laurent polynomials [T^ EU]. 
Moreover, as pointed out by Kharchev et al. pTj and Suris [22], the Ablowitz- 
Ladik hierarchy is equivalent to the relativistic Toda hierarchy. It is rather 
the relativistic Toda hierarchy that is directly connected with the 2D Toda 
hierarchy. Brini et al. [23j characterized this integrable hierarchy (which 
they call the Ablowitz-Ladik hierarchy) as a kind of "rational reduction" of 
the 2D Toda hierarchy. The notion of rational reduction was first introduced 
in the case of the KP hierarchy [2l] . In the Toda version, the Lax operators 
L and L of the reduced system are expressed as "quotients" of finite order 
difference operators in the variable s. 

We show that the Lax operators of the solution in question do have such 
a factorized form. To this end, we consider a factorization problem of Z x Z 
matrices [251 ES EZ] that can capture all solutions of the 2D Toda hierarchy. 
In general, solving this factorization problem directly is extremely difficult. 
In the present case, however, we can find an explicit form of the solution at 
the "initial time" t = i = 0. This is enough to determine the "initial values" 
of the Lax operators. As expected, the initial values turn out to have the 
factorized form of Brini et al. Since the factorized form is preserved by the 
flows of the 2D Toda hierarchy (this is what the notion of reduction means), 
we can conclude that the solution in question belongs to the Ablowitz-Ladik 
or relativistic Toda hierarchy. 

Let us emphasize that shift symmetries are interestesting in themselves. 

"'^The bars " ~ " do not mean complex conjugation. 
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This is another subject of this paper. We shall encounter a few new aspects 
in the course of the consideration on the integrable structure. Firstly, we 
derive the shift symmetries more carefully than in our previous work [H |2] . 
This explains why a careless use of these symmetries leads to a contradictory 
result |5]. Secondly, we translate the shift symmetries to the language of 
Z X Z matrices. Remarkably, this reveals a close relationship with quantum 
dilogarithmic functions j2Sl 122] • In the matrix representation, the vertex 
operators [9l [H] in the fermionic expression of Z{s,t) and Z'[s,t,i) turn 
into matrix-valued quantum dilogarithmic functions. The shift symmetries 
thereby become a straightforward consequence of an infinite product formula 
of the quantum dilogarithmic functions. We believe that this issue deserve 
to be studied in a more genera context. 

This paper is organized as follows. Sections 2, 3, and 4 presents the per- 
spective from fermions and tau functions. In Section 2, the partition func- 
tions of the two melting crystal models are formulated. The setup of the 2D 
complex free fermion system is also explained here. In Section 3, the quantum 
torus algebra of fermion bilinears is introduced, and the extended set of shift 
symmetries are derived. These algebraic relations are used in Section 4 to 
derive the tau functions from the deformed partition functions of the melting 
crystal models. Sections 5 and 6 are mostly concerned with infinite matrices 
and Lax equations. In Section 5, the shift symmetries are reformulated in the 
language of infinite matrices. This is a place where matrix-valued quantum 
dilogarithmic functions show up. Although the matrix versions of shift sym- 
metries are not necessary for the subsequent consideration, the proof of these 
algebraic relations overlaps with technical details of the contents of Section 6. 
Section 6 is devoted to the proof of our conjecture that the modified melting 
crystal model gives a solution of the Ablowitz-Ladik hierarchy. The relevant 
matrix factorization problem is formulated and solved here. The associated 
Lax operators are thereby shown to take a factorized form as expected. 

2 Melting crystal models 
2.1 Previous model 

The partition function of the undeformed model is the sum 



over the set V of all partitions A = (Ai, A2, • • •)• q and Q are complex numbers 
with IQI < 1 (or may be thought of as formal variables). sx{q~^) is the 




(2.1) 
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special value of the Schur function sx{x) of infinite variables x = {xi,X2, ■ ■ ■) 
at 

This special value has the well known hook formula [30] : 

where k(A) are the second Casimir value (equivalently, twice the total con- 
tent) 

oo 

'^(A) = 5Za,(A,-2^ + 1) = 2 U-') 
i=i (i,i)eA 

and h{i,j) is the hook length of the cell in the Young diagram of shape 
A. By the Cauchy identity, 

sa(xi, X2, . . .)sx{yu 2/2, • • •) = H " ^i2/i)~\ (2-3) 

\€V i,j>l 

one can rewrite the sum Z to an infinite product: 

oo oo 

^ = n (1 - Q^'^'-'r' = 11(1 - Qg")-". (2.4) 

i,j=l n=l 

We deform (12.1 p with a discrete variable s G Z and an infinite number of 
continuous variables t = (ti,t2, ■ ■ ■) as 

Z{s, t) = Y, SA(g-0'Q'^'+'^'^'^^'e*(^'^'*), (2.5) 
xev 

where $(A, s,t) is a linear combination 

oo 

$(A,s,t) = ^tfc$fc(A,s) 

k=l 

of the external potentials 

oo oo 

$^(A, s) = ^ ^fc(A,+s-i+i) _ ^ (formal expression) 

i=l i=l 

oo 1 fcs 

= ^(g^(^»+^-^+i) - g^(-^+i)) + ^^^g^ (2.6) 
1=1 ^ 

The formal expression in the first line is related to the normal ordering pre- 
scription in the fermionic expression (I2.14p of these potentials. 
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2.2 Modified model 



In the modified model, one of tlie factors of the weight is modified as s\{q^^) — )• 
stx{q~P), where *A is the conjugate (or transposed) partition of A. Thus the 
undeformed partition function reads 

By the dual Cauchy identity 

oo 

sx{xi, X2, . . .)s txivu 1/2, • • •) = n (^ + ^^^j)' (2-8) 
xev i,j=i 

Z' can be cast into an infinite product: 

oo oo 

Z'=1[{1 + Qq'^^-') = + Qq^r- (2-9) 

i,j=l n=l 

Since stx{q^^) and sx{q~'') are related as 

s^xiq-n = q^^^^^'sxiq~n- 
by the hook formula (12.21) . one can rewrite the sum (12.70 as 

Z' = ^SA(g-OV^^^^'Q'^'- (2.10) 
Xev 

Thus Z and Z' are the first two members of the infinite family 

^i = $^SA(g-OV"^'^/'Q'", / = o,i,..., (2.11) 

Xev 

of generating functions. Zi is a generating function of local Gromov-Witten 
invariants in the local CP^ geometry Xi = 0{-l) © 0{l - 2) ^ CP^ ^Tj. Z 
and Z' can also be derived by the method of topological vertex ^ as the 
amplitudes of topological string theory on Xq and Xi. Note, however, that 
the parameter Q in the string amplitude of Xi is flipped to —Q [8J. This 
difference affects the final result (see Section 6) on the underlying integrable 
structure. 

We now introduce two sets of continuous variables t = {ti,t2, ■ ■ ■) and 
i = {ii, ^2, • • •)) deform the sum (12. 7p as 

Z'is,t,i) = ^SA(g-05u(g-0Q'"+'^'+'^/'e*(^'''*'*\ (2.12) 
Xev 
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where $(A, s, t, t) is the hnear combination 

oo oo 

$(A, s, t,t) = Y^ tfc$fc(A, s) + 4$-fe(A, s). 
fc=i fe=i 

A)'s are defined by the same formula as (12.61) with k replaced by —k. 
2.3 Fermionic formulation 

We use the same formulation of 2D complex free fermions as in the previous 
work [HE]- The Fourier modes 'i/'n, V'*, n G Z of the fermion fields 

satisfy the anti-commutation relations 

The Fock spaces are decomposed to charge-s sectors {s G Z) with the ground 
states 

(s| = (-00| ■ ■ ■ tfj^-l^t^ \s) = ^-si^-s+l ■ ■ ■ I - oo) 

and the excited states 

(S, A| = (-00| ■ ■ ■ ^^2+.-l^Ai+s' = ^-Ai-.V^-A2-s+l • • ■ I - oo) 

labelled by partitions. 

Let us introduce the fermion bilinears 

Lo = ^n:tlj_n^l:, Wo = ^n^-.i/j-n^l:, 

and the vertex operators P, [H] 

r±(z) = exp I ^ — J±fc J , r'^(2) = exp I - ^ — ^ J±fc 
\fc=i / \ fc=i 

It is convenient to extend the definition of the vertex operators to multi- 
variables a; = (xi, X2, . . .): 

r±(a;) = nr±(x.), r^(a.) = nr'±(x.). 

i>\ i>\ 
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The matrix elements of these multi- variable vertex operators with respect to 
the states (A,r|, \fi,s) are skew Schur functions [12j: 

{X,r\r_{x)\fi,s) = {fi,s\r+{x)\X,r) = 6rsSx/f,{x), 
{X,r\T'_{x)\n,s) = {n,s\T'^{x)\\,r) = 5rsS tx/ t^,ix). 

Thus the special values of the Schur functions in the definition (12.51) and 
fl2.12p of the deformed partition functions can be expressed as matrix ele- 
ments of the vertex operators specialized to a:; = q^''. Since the potentials in 
the partition functions, too, can be expressed as 

\X\ + s{s + l) = {X,s\Lo\X,s), MX,s) = {X,s\Hk\X,s) (2.14) 

and the other matrix elements of Lq and Hk vanish, one can rewrite the 
partition functions as 



and 



where 



Z{s,t) = (.|r+(g-Og^°e^Wr_(g-'')|s) (2.15) 



Z'{s,t,i) = (.|r+(g-OQ^"e^(*'*)r (g-^k), (2.16) 



H(t) = ''^tkHk, H(t, t) = tfcfffc + tkH_k. 

k=l k=l k=l 

3 Shift symmetries in fermionic formulation 

3.1 Three sets of shift symmetries 

Let us introduce another set of fermion bilinears {Vm ''}k,mez'- 

J\z\=R 27rz f-^ 
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The radius i? > is chosen arbitrarily. Hk and Jk are particular cases of 
these fermion bilinears: 

Vm s satisfy the commutation relations 

[Vi')^^)] = (,a-M/2 _ ^(.n-M/2)(V^j^+0 _ 5^^„,o^^) (3.1) 



for k and I with k + I and 

[Vl'\ Ki"')] = - + rn5^+n,o, (3.2) 

thus give a reahzation of (a Central extension of) the 2D quantum torus Lie 
algebra 

"Shift symmetries" connect these fermion bilinears. There are three sets 
of fundamental shift symmetries: 

Theorem 1. 

(i) For k > and m E Z, 



(3.3) 



(ii) For k > and m G Z, 



(3.4) 



(iii) For k,m E 1j, 

f l3.3p and fl3.5p are presented in our previous work [H |2] . fl3.4p is new, and 
we show the derivation of these shift symmetries in detail below. We shall 
also see why (13. 3p and (13. 4p do not hold for A; < 0. 

Remark 1 . In our previous work [H [2] , shift symmetries are presented in such 
a form as 

(3-6) 



-1)" ( Ki+fc ~ l-qk ^'"+^'0 



and 



, ^ ^ (3.7) 



In this paper, however, we dare to avoid (13.61) and fl3.7l) because, firstly, it is 
fl3.3p and f l3.4p that are used in the subsequent consideration and, secondly, 
the operator products r_(g-^)r+(g-^) and T'_{q-P)T'^{q'P) in and ([32]) 
require a careful treatment. 



3.2 Derivation of shift symmetries 

Derivation of (13. 4p is parallel to the proof of (13. 3 p [H [2] . We start from the 
equality 

■.ip(z)ilj*(w): = ip(z)i/j* (w) 

z — w 

that holds for \z\ > \w\. Since |g| < 1, this equality can be specialized to 

k/2 



q" 



{l-q'')z' 

The fermion bilinears Vm in (13. 4p can be thereby expressed as 

VL-'^ = q-'/' i ^.z-^lj{q-^"z)r{q"'z) - -^5^,0 
J\z\=R ^TTi 1 - g'^ 

We now use the formulae 

T'^{x)ij{z)T'^{x)-^ = (1 + 
T'^{x)i)*{z)T'^{x)-^ = (1 + X2^^)-V*(^) 



(3.9) 



that holds when < 1. The formulae for T'_^{x) in (13. 9p imply that the 

equality 

rV iq-''Mq-''^^z)^* (g^/'2)rV (q-P)-^ 

k 
i=l 

holds for \z\ < Multiplying this equality by g~^/^z'"/27rz and inte- 

grating it along the circle \z\ = Ri < Igl*^^"^-*/^, we find that 



Xq-T' 

r dz ^ ^^-^^^ 
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On the other hand, the formulae for T'_(x) in fl3.9p imply that the equality 

T'_{q-P)-^^{q-''^h)^*{q''/h)r_{q"P) 



k 



i=l 



holds for \z\ > |gp -"^^Z^. Note here that the prefactors of tp{q ^/'^z)il)*{q^/'^z) 
in this equality and the previous ones are connected by the simple relation 



fc k 

m 



Z 

i=l 



Y[{1 + q'-O'+^y^z) = J](l + (3.11) 



Therefore, multiplying this relation by q ^/'^z'^^^ and integrating it along 
the circle \z\ = R2 > Igp"^-*/^, we obtain the equality 



r-iq-T' f vl;? + ^W,o ) r (g- 



-k/2 



J\z\=R, 2m A A 



(3.12) 



The right hand side of fl3.10p and 03.121) are almost the same except for 
the difference of the contours. Actually they coincide, because the integrands 
have no singularity in the annulus Ri < \z\ < R2, and the contour can be 
continuously deformed therein without changing the contour integral. This 
completes the proof of f l3.4p . 

3.3 Where shift symmetries break down 

Let us consider what occurs when < 0, in other words, when k is replaced 
by —k in the setting of the proof of f l3.4p . We start from the equality 

:ilj{z)ip*{w): = —'iIj*{w)iIj{z) 



w — z 

r{k) 



that holds for \w\ > \z\. The fermion bilinear Vm can be thereby expressed 
as 

Vi'^ = ~q"' / ^z^r{q-^'^z)i;{q'/'z) + -^5^,0- (3.13) 
J\z\=R 27^^ 1 - q'' 
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Under the adjoint action of r'_j_(g ''), the operator product ip*{q ^/'^z)il){q^/'^z) 
transforms as 

T'^{q-nr{q-'^'z)i,{q'/'z)T\{q-'')-' 
= ri(l + q'-^^'^^^'^z)-^ ■ i,\q-''l''z)i,{q^l\) 



i=l 



and 

r_{q-P)-^^*{q-^/^z)^lj{q^/^z)Ti{q-f) 

k 
i=l 

Consequently, we obtain the following analogues of fl3.10p and fl3.10p : 

^'-.(q-n [-yJn^ + r^k^rnfi^ Kiq-T' 

= q'^' / ^z"^ 17(1 + ■ ^P*{q~'/'z)^P{q'/'z), 



Unlike the proof of (13 ■4p . the integrands of these contour integrals have 
poles, one of which sits on the circle 1^1 = \q\^^~^^/'^. The contours \z\ = Ri 
and \z\ = R2 are separated by this circle, and cannot be reached from the 
other by continuous deformation without crossing the poles. Thus an equality 
like (13 ■4p does not hold. 

This explains why (13. 4p breaks down for k < 0. By the same reasoning, 
one can see that fl3.3D does not hold for < 0. 
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4 Partition functions as tau functions 
4.1 Previous model 

In our previous work [U |2], the partition function Z{s,t) is shown to be 
related to a tau function t{x, t) of the ID Toda hierarchy as 



Z{S, t) = exp ( ^ ) g-(^+l)(2^+l)/6^(,^ (4.1) 



where L{t) denotes the alternating inversion 

i{t) = {-h,t2,-h,..., {-!)%,...) 
of t = (ti, t2, . . .). The tau function t{s, t) has the fermionic expression 

T{s,t) = {s\exp ^^^^tkJk^ g\s), (4.2) 

where 

g = g^o/2r_(g-or+(g-OQ''°r_(g-or+(g-Og'^°/'. (4.3) 

The two sets (13. 3p and (13. 5p of shift symmetries are used to derive (14. ip from 
the fermionic expression (12.151) of the partition function. 

By the same same shift symmetries, one can also show that g satisfies the 
algebraic relations 

Jkg = gJ-k, A; = 1,2,.... (4.4) 
This implies that the associated tau function [13] 



r(s,t,t) = (s|exp l^^^tkJkj gexp y-^^ttJ-kj \s) 

of the 2D Toda hierarchy depends on t and i through their difference: 

r(s, t, i) = r(s, t — i). 
Consequently, r(s,t) can be expressed in different forms, e.g.. 



r(s,t) = (sl^fexp I ^tkJ^k I \s) 

(4.5) 



vfc=l 



[s \ exp 



^^tkJkjg [\^tkJ-k] \s). 
, fc=i / \ k=i ) 
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4.2 Modified model 

We now turn to the modified model, and derive an analogue of (14.11) from 
the fermionic expression fl2.16p . The derivation is mostly parallel to the case 
of the previous model [H |2]. 

First step: The first set (13.31) of shift symmetries, applied to Hk = Vq''\ 
k = 1,2, . . ., imply that 

Moreover, the third set (13. 5p of shift symmetries imply that 

yi^) = q-Wo/2j^^Wo/2_ 

From these equalities, we obtain the algebraic relation 

k 

T+iq-'')HkT4q-P)-' = i-l)'T4q-^r\-'^^/'Jkq'^<^/'T4q-^) 



q- 



l-q^ 

that connects and J^. This algebraic relation can be exponentiated as 



q^t 



V^{q-P) exp Y^^uHu T^^q-^)-' = ^xp ^ 



k 



.fc=l / \A:=1 ^ 



X 

(4-6) 

X r_(g-'')-ig-^°/2exp [YS-l)%Jk\ q'^'/^T _{q-P). 

Second step: The second set (13. 4p of shift symmetries, applied to H_k = 
Vq^''\ k = 1,2, . . ., imply that 

T'_iq-n-'H_,T'_iq-P) + ^ = T'4q-'^)v[l'^T'4q-T'- 

Moreover, the third set (13. 5p of shift symmetries imply that 

yi-k) ^ ^-Wo/2j_^^Wo/2_ 

Consequently, we obtain the algebraic relation 



h ' 
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which can be exponentiated as 

T'_{q-T'exp (f^hH^ r = exp (-f^^L_ 

\k=i J \ fc=i ^ 



X 

(4.7) 



Third step: Let us factorize the operator in fl2.16p as 

T+{q-'')QV^'~'^r_{q-'') = r+(g-Oexp (^fj^'^^'^j 

X r+(g-OQ^"r'_(g-Or'_(g-0-'exp (f^hH^k] ^ ^0- 
By (14. 6p and fl4.7l) . we can rewrite it as 



X 



vfc=i 



q ^°/2x 



.fc=i 

oo 



X 



where 

g' = g^o/2r_(^-P)r^(^-P)gLor'_(^-P)r;(g-p)g-^o/2_ (4,^ 

Plugging this expression into (12.161) and noting that 

we end up with the following result: 
Theorem 2. Z'{s,t,i) can be expressed as 



oo 



Z'{s, t, i) = exp E 'i^^ r'{s, .(t), -i), (4.9) 
where T'{s,t,i) is the tau function 

T'{s,t,t) = (s|exp I ^tkJk j fi-'exp | - ^hJ-k ) |s) (4.10) 

\fc=l / \ k=l 

of the 2D Toda hierarchy defined by the operator 
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Remark 2. If the shift symmetries fl3.3p and (13 ■4p were vahd for k < as 
well, one could show that g' satisfy the algebraic relations 

J±kg' = g'J±k, k = 1,2, . . . {wrong statement) . 

This imply that the tau function would become an almost trivial one: 

T'{s,t,t) = exp I — "^^ktktk j {s\g'\s) {wrong statement). 

\ k=l J 

This is a typical example of wrong consequences derived by a careless use of 
shift symmetries [Sj. 



5 Shift symmetries in matrix formulation 

5.1 Matrix representation of fermion bilinears 

It is well known [12J that the mapping 

X = ^ XijEij t-> X = ^ Xifip-iip*: (5.1) 

gives a projective representation (in other words, a central extension) of the 
Lie algebra gl(c>o) of Z x Z matrices. In the following, we naively identify 
X and X. The fundamental fermion bilinears Lq, Wq, H^, can be thereby 
identified with the Z x Z matrices 

Lo = A, W, = /\\ Hk = q''^, Jfc = A^ (5.2) 

where A and A are the matrices 

A = ^ iEii, A = ^ Ei^i+m 

that satisfy the twisted canonical commutation relation 

[A,A] = A. (5.3) 

A and A amount to the shift and multiplication operators e^l^^ and s in the 
Lax formalism of the 2D Toda hierarchy [TOl [TT| . 

The fermionic realization {Vm ^}k,m& of the quantum torus Lie algebra 
correspond to the Z x Z matrices 
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These matrices form an associative algebra rather than a Lie algebra. The 
commutation relations 

can be derived from the fundamental relation 

Ag'^ = gg^A (5.6) 

of the generators A and q^. 

Last but not least, the vertex operators r-i-(x) and r'_j_(x) turn out to have 
a matrix representation of the form 

00 u 
X 



T^{x) = exp I ^ -k^' ) = (1 - xA±i)- 



00 / \ I, 



(5.7) 

r'±(x) = exp [ - ^ ^^^^ ) = (1 + a^A^')- 

The inverse of 1 — xA^^ is understood to be the geometric series 

(1 - xK^^)-' = 1 + xA±i + x^K^^ + ■■■ . 

Among the four matrices of (15. 7p . those with label "+" are upper triangular, 
and those with label "— " are lower triangular. The matrix representation of 
^±{<1~^) and r'_j_(g~'') is an infinite product of these matrices specialized to 

X = 

oo 



Thus r-|-(g P) and r'_j_(g p) may be thought of as matrix- valued quantum 
dilogarithmic functions in the sense of Faddeev et al. [28l [29] . 

5.2 Reformulation of shift symmetries 

The shift symmetries (13. 3p . (13. 4p and (13. 5p can be reformulated in the matrix 
representation as follows. 
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Theorem 3. The Z x Z matrices / 1 5.^)) . ( [5.^[ j anc? ( (5. gj) satisfy the following 
analogues of / f5^) . anc? / f5^) .- 

r+r'')VWr+(g-0-^ = {-lfT^{q-r'vi::lT_{q~n, (5.9) 

iiTere A; anc? m in \5. 9\) and li5.1U\) range over k > and m G Z, and ( I5.ii)) 
/ioWs /or k eTj and m G Z. 

Proof. For a comparison with the proof of (13.41) . let us show the derivation 
of flS.lOp in detail. A technical clue is the algebraic relation 

g-'^^AV^ = /"A" (5.12) 

that is a partially exponentiated form of (15. 3p . 

We first consider the left hand side of (I5.10p . namely, 

(I5.12P implies that 

oo oo 
j=l i=l 

Multiplying this equality by r'^(g~^) and g"'^^ from the left and right sides 
yields the equality 

oo oo 

r'_,(g-Og"'^r;(g"0"' = + ?'"'^'A) ■ J](l + g^+'^-i/^A)-! ■ g-'^^ 

i=l i=l 
k 

= n(i + g^-'/'A)-g"'^, 
j=i 

which implies that 

k 
i=l 

We now consider the left hand side of (I5.10p . namely, 

^'-{q'T'yi'+k^'-iQ'") = g'=(-+^)/2A™+fcr'_(g-0"'g"'^r'_(g-0- 
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In the same way as we have followed above, we can derive the equality 

oo oo 

i=l 1=1 
k 



i=l 

Since 



^fcV2Afc ^^(i + g^-fc-i/2A-i) = + q'-y^A), (5.13) 



i=l i=l 

we can conclude that 



T'_{q~T'vi-_,'iT'_{q~P) = ^^2^- + q^-^^A) . q 



kA 
i=l 



Thus both hand sides of (15.101) turn out to coincide. We can derive (I5.9p 
in the same way from the equalities 



^+{q'')Vi^^T+{q-P)-' = g-^™/2A™ JJ(1 - q'-'^+^^A) ■ g^^, 

i=l 

k 

r_(g-0-Vj2,r_(g-0 = q~Kn^+k)/2^m+k "Q^^ _ g^-i/2A-i) . q 



i=l 



and 

k 



'A^ JJ(1 - q'-y^A-') = (-l)'=g-W2A™ JJ(i _ g^-fc+i/2A) . 



A;=l fc=l 



Lastly, (15.111) can be derived by straightforward calculation of the matrix 
elements of both hand sides as done in the language of fermions [HE]- □ 

5.3 Where shift symmetries break down 

Just as in the fermionic formulation, (15. 9 p and (I5.10p do not hold for k < 0. 
Let us replace k ^ —k in (15. 9p and (I5.10p . and show that the equalities do 
not hold. 

By replacing k — )■ —k, both hand sides of (15. 9 p turn into 

T'{q-P)Vi^^T'4q-P)-^ = q-''"'^^A'^T'{q-P)q^^T'4q-f')-^ 
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and 

T'_{q-Py^V^lf,T'4q-'') = q-''^'''-''^/^A"'-''T'_{q-P)-'q''^T'_{q-P). 

As in the proof of Theorem 3, we can calculate the triple product on the 
right hand side as 

k 

Ki<l~')Vi^^^'+iq~')'^ = q-'^'^I^K^ + q'-'^-^'^KY^ ■ q^^ 

i=l 



and 

k 

i=l 

Let us consider the inverse matrices (l+g*-'^-^/^^)-^ and (l+g*-^/^^^^)-^ 
that show up in these expressions. They are understood to be given by 
geometric series: 

(1 + q'~'^^A-y = 1 - q'-'/^A-' + iq'-'/^A-')^ 

In particular, (1 + ^^-^-1/2^)-! and (1 + q'-^/'^A"^)-'^ are upper triangular 
and lower triangular, respectively. Consequently, unlike the equality fl5.13p 
that underlies flS.lOp . we have the inequality 

k k 
i=l i=l 

fl5.14p implies that 

T'4q-nVi'^T'4q-r' ^ r'_{q~r'vi%T'_{q-n. 
By the same reasoning, one can confirm that 

TAQ-')vi-'^T4q~r' ^ (-i)'r-r'')-virir_(g-o. 

6 Relation to Ablowitz-Ladik hierarchy 

6.1 Two reductions of 2D Toda hierarchy 

In the Lax formalism [TOl [TT], the 2D Toda hierarchy is formulated as a 
system of Lax equations 

^ = [Bk,L], |^ = [5fc,L], 
otk oik 
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for two Lax operators L and L. In the case of a bi-infinite lattice, L and L 
are difference operators (i.e., linear combinations of the shift operators e"^'' 
in the lattice coordinate s) of the form 



n=l n=l 



where m„ and Un depend on s, t and t. and B^ are constructed from L 
and L as 

where ( )>o and ( )<o denote the projections to linear combinations of e^^" 
for n > and n < 0. 

The ID Toda hierarchy and the Ablowitz-Ladik hierarchy are two typical 
reductions of the 2D Toda hierarchy. Let us recall the reduction procedure. 

Reduction to ID Toda hierarchy 

The ID Toda hierarchy can be derived from the 2D Toda hierarchy by im- 
posing the reduction condition 

L = L-^. (6.2) 
Both hand side of this equality define a difference operator of the form 

£ = + 6 + ce-^' (6.3) 
that satisfies the Lax equations 

Actually, since Bk, Bk and £ are related as 

Bk + Bk — 

the time evolutions in t and i are no longer independent in the sense that 
d£ 02 ^ 

;S + sC = '''"'="^'^"'=" = °- 

Thus the time evolutions in i are redundant, and we are left with just one 
set of Lax equations 

= S, = (£^)>o. (6.4) 

Otk 

This is one of the well known Lax forms of the ID Toda hierarchy. 



21 



Reduction to Ablowitz-Ladik hierarchy 

The Ablowitz-Ladik hierarchy can be derived from the 2D Toda hierarchy as 
a kind of "rational reduction". In the formulation of Brini et al. [23], L and 
L are chosen to be quotients 

L = BC-\ L-^ = CB-^ (6.5) 

of two difference operators of the form 

B = e^'-b, C=l-ce~^\ 

B~^ and are understood to be difference operators of the form 



B-' = - J2ib''e^')''b'^ = -b{s)-^ - J2 K^y^ ■■■b{s + k)-^e 

k=0 k=l 
oo oo 

C-^ = l + J2ice~^'f = l + J2 ^(^)^(^ - 1) ■ • ■ c(s - + l)e-' 



(6.6) 



k=l k=l 



where b{s) and c(s) are abbreviations of c(s, t, t) and c{s, t,t). B and C have 
inverses of different types, 



B-^ = J2 e"^' (be-^^)'' = 6"^= + ^ 6(s - 1) ■ ■ ■ 6(s - A; + 1) 

fc=0 k=l 

oo oo 

= - ^(e^=c-^)^ = - ^ c(s + 1)-^ • ■ ■ c(s + k)-^e 



k=l k=l 



but they are not used here. (16.61) imply that (CB^^)^'^ ^ BC^^ . Thus we 
can avoid the trivial case where L = L = e^" . Under the reduction condition 
(16.51) . the Lax equations (16. ip can be reduced to the following equations: 



Qf^ = {{BC-'f)>,B-B{{C-^Bf)^_^ 

^={{BC-r),,C-C{iC-'Bf)^^ 
^'^ (6.7) 
{{CB-')')^,B-B{{B-'C)'), 



dB 



^^ = {{CB-'f),,C-C{{B-'Cf)^^. 
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These equations have another expression of the following form: 



{{BC- 




H B {{C- 




{{BC- 


")')<o^- 


fC((C- 




{{CB- 




^B{iB- 




{{CB- 




-C{{B~ 





dB _ 
dtk 

_ 

dtk 
dB _ 

dtk 

_ 

dtk 

We can deduce from these two expressions of (16.71) that the right hand side 
are difference operators of the form fk, gke~^% fk and gkC'^" . Thus (16. 7p can 
be reduced to evolution equations of the form 

db _ dc _ 9b _ - dc __ 
dtk dtk dtk dtk 

This implies the consistency of reduction. In particular, the factorized form 
(16. 5p of the Lax operators is ensured to be preserved by the flows of the 2D 
Toda hierarchy. (16. 8p can be further converted to the usual Hamiltonian form 
[22] of the Ablowitz-Ladik hierarchy by a change of variables [23] . 

Another formulation of reduction to Ablowitz-Ladik hierarchy 

The Lax operators of (16. 5 p have another expression 

L = C-^B, L-^ = B-^C, (6.9) 
where B and C are difference operators of the form 

B = e^^ -b, C=l- ce-^' 
that satisfies the relation 

CB = BC. (6.10) 

B~^ and C*~^ are understood to be defined by series expansions hke (16.61) . 
The operator relation (I6.10p turns into the functional relations 

b{s) + c{s) =b{s) + c{s + l), c{s)b{s-l) = b{s)c{s). (6.11) 
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b, c and b, c are thus mutually connected by an invertible transformation. 
The Lax equations (16.1 p can be reduced to the equations 

dB (,~~ 



otk V / >o V / >o 



dC 



at. ((^'^")')>o^-^(('^"^''y>„ 



dB (~~ 



(6.12) 



'-^^^[m-^f)j-c{^B'^cf)^^ 

for B and C. These equations can be further converted to evolution equations 
for b and c. 

The factorized form (16. 9p is suited for a comparison with the relativistic 
Toda hierarchy [IS]. If the Lax operators are factorized in this form, the 
eigenvalue problem 

m = z^, L-^^ = z~^^ (6.13) 



for the Baker- Akhiezer functions ^ and of the 2D Toda hierarchy ITTj 
can be converted to a "generalized eigenvalue problem" of the form 

m = zC^, M = zC^. (6.14) 

As pointed out by Kharchev et al. ^J, this generalized eigenvalue prob- 
lem can be derived from the traditional 2x2 matrix Lax formalism of the 
Ablowitz-Ladik hierarchy as well. One can formulate an auxiliary linear 
problem of the relativistic Toda hierarchy on the basis of this generalized 
eigenvalue problem. 



6.2 Factorization problem of infinite matrices 

Difference operators in the Lax formalism of the 2D Toda hierarchy are in 
one-to-one correspondence with Z x Z matrices by the rule 

oo oo 

J2 «fc(s)e^^' ^ Yl diag[afc(s)]A^ 

fc=— oo k=—oo 

where diag[a(s)] stands for the diagonal matrix 

diag[a(s)] = ^ a{i)Eii. 
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Since this correspondence preserves the algebraic structure, one can reformu- 
late the Lax formalism of the 2D Toda hierarchy in terms of Z x Z matrices. 
Let us use the same notations L, L, 5„ etc. for the corresponding Z x Z 
matrices. 

In this matrix formulation, one can formulate a matrix factorization prob- 
lem that can capture all solutions of the 2D Toda hierarchy [2^1 EZ] • Given 
a constant invertible Z x Z matrix U, the problem is to find two Z x Z matrices 
W = W{t,t) and W = W{t,t) with the following properties: 

(i) is a lower triangular matrix, and all diagonal elements are equal to 
1. 

(ii) W is an upper triangular matrix, and all diagonal elements are non- 
zero. 

(iii) W and W satisfy the factorization relation 

exp ^XJ^feA^^j ^^^^ (^"XJ^'^^^'^j = W-^W. (6.15) 
If such matrices W, W exist, they satisfy the so called "Sato equations" 

(6.16) 



dW I. dW - 

-- = BkW -WA.\ =BuW, 
dtk otk 



dW - dW - - - I. 
-- = B,W, =B,W-WA-\ 
otk otk 

where Bk and Bk are constructed from W and W as 

Bk = [WA'^W-') , Bk= {WA~'W-') . 

( )>o and ( )<o stand for the upper triangular part and the strictly lower 
triangular part. The Sato equations 06.161) imply that the matrices L and L 
defined by 

L = WAW-\ L = WA-^W-^ (6.17) 

satisfy the Lax equations (16. ip . 

The factorization problem (16.151) is directly connected with the fermionic 
expression [T3] 



oo 



T{s,t,t) = (s|exp ^tkJk fi'exp -^tkJ-k \s) 



vfc=l / \ k=l 
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of tau functions. Namely, the Z x Z matrix U = {uij) corresponds to the 
operator g by the Bogohubov transformation 

9^j9~^ = ^'ipiUij (6.18) 

of fermion creation-annihilation operators [12]. Its infinitesimal form is ex- 
actly the correspondence (15.11) between Z x Z matrices and fermion bilinears. 

The operators (14. 3 p and (14.81) defining the tau functions r(s, t) and r'(s, t, i) 
correspond to the matrices 

U = g^'/2r_(g-0r+(g-0Q^r_(g-0r+(g-0g^'/' (6.19) 

and 

u' = g^'/2r_(g-or+(g-'')g^r'_(g-orV(g-'')g-^'/', (6.20) 

where r±(a;) and r'_j_(a;) are the matrix representation (15.81) of the multi- 
variable vertex operators. We now consider the factorization problem (I6.15p 
in these particular cases. 

6.3 Initial values of Lax operators 

When t and i are specialized to t = t = 0, we can find an explicit form of 
the solution of the factorization problem (I6.15P for the matrix (16.200 : 

Lemma 1. The factorization problem 

U' = W-^W (6.21) 
for the matrix 116. 20\} can be solved as 

W = g^'/2r'_(gg-p)-ir_(g-^)-ig-AV2^ 

Proof. We can use the algebraic relations 
to rewrite the triple product in the middle of (I6.20p as 

r+(g-OQ^r'_(g'0 = Q^r^{Qq~^)r'_{q-n 

= Q^Ti{q-P)T^{Qq-n 
= T'_{Qq-nQ^T+{Qq-P). 



(6.22) 
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Thus U' becomes a product of triangular matrices: 

We have only to modify the first (lower triangular) part with the diagonal 
matrix q~^ so that all diagonal elements become 1. □ 

The initial values of the associated Lax operators f l6.17p at t = t = are 
obtained from (16.221) as 

L = WAW~\ L-^ = WA'^W-\ 

These matrices turn out to be a quotient of two matrices: 

Lemma 2. 

L = (A - g^)(l + Qq^-'A-y = (1 + Qq^A-y\A - q^). (6.23) 
Proof. By construction, L is a matrix of the form 

L = q^'/'TLiQq-T'^-iq-T'q'^'^'^q^'^'TAq'n^'-iQq'nq~^'^'- 

Let us calculate this matrix step-by-step. Note that specialization of (15. lip 
to /c = reads 

This implies that 

The next stage is to calculate the product 

r4Qq-p)-'^4q-p)-'q'^'^'^q^'^^^-iq~n^'-iQq'n 

= AT'_ {Qq-T'^- {q-n-\^-'/^T_ {q-^)Ti (Qq-f) 

with the aid of (I5.12p . Firstly, the triple product of r_(g~'')^^ and q^ can 
be calculated as follows: 

r_(g-0-^g^r_(g-0 = T4q-n-' ■ q^T^q-^q'^ ■ q^ 

oo oo 

= - g-i/^A-i) . J](l - g^+i/2A-i)-i . 

i=l i=l 

= (1 - q'/^A^')q\ 
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Since 1 — g^/^A ^ commutes with T'__{Qq ^)^^, the next task is to calculate 
a triple product of T'^i^Qq"^)"^^ and q^. This can be achieved as follows: 

T'_{Qq-n-'q^r'_iQq-^) = T'_iQq-^)-' ■ gV^Qg-^g"^ ■ g^ 

oo oo 

= + Qq'-'/'A-y ■ + Qq'+'/'A~') ■ q^ 

i=l i=l 

= (i + ggV2A-i)-v_ 

These calculations yield the equality 

r'_ {Qq-T^T_ {q-P)-\-^'I^Xq^'I^T^ {q~'')T'_ {Qq-") 
= A(l - gi/2A-i)(l + QgV2A-i)-igA-i/2^ 

hence an expression of L of the form 

L = g^'/2(A - gV2)(l + ggl/2^-l)^A-l/2g-AV2_ 

We now use the algebraic relations 

^AV2^^-AV2 ^ g-A-1/2^^ ^AV2^-1^-AV2 ^ 

Ag^ = g^+^A, A'^g^ = g^-^A^^ 
that follow from fl6.24l) and (15.121) to rewrite L further as follows: 

L = (g-^-V2A - gV2)(l + g^A^-l)^A-l/2 

= g^-^/'(g-^+^/'A - gV2)(i + QgA-iA-i)-i 
= (A-g^)(l + Qg^-iA-i)-i. 

This gives the first part of (16.231) . In the same way, starting from 

L = q^'/\l + Qq^/^K-')-\k - q'/')q^-'/\~^'/' 

and repeating similar calculations, one can derive the second part of (16.231) . 

□ 

Lemma 3. 

I-i = (1 + gg^-iA-i)(g^ - A)-i = (g^ - A)-\l + Qg^A-^). (6.25) 

Proof. The proof is parallel to the case of the previous lemma. By construc- 
tion, is a matrix of the form 

= g^^/2Q^r+(gg-0rV(g-0g-^'/'A-ig^'/2 

X T'^{q-T'rAQq-T'Q-''q-^'^'- 
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With the aid of fl6.24p and fl5.12p . one can derive the equahty 

This imphes that 

Z-i = q^'/'Q^{l + g-V2A-i)(i _ g^i/2^)-i^-Ag-A^-AV2_ 

It is easy to derive the final form fl6.25p of from this expression. □ 
6.4 Identification of solution 

Translated to the language of difference operators, the foregoing results f l6.23p 
and fl6.25p show that the initial values of the Lax operators at t = t = are 
factorized as 

L = BC-^ = C-^B, L-^ = -CB-^ = -B-^C (6.26) 

with factors of the special form 

B = B = e^^ -q\ C = 1 + Qq'~^e'^^ , C = 1 + Qq'e~^^ . (6.27) 

Note that this expression is slightly different from (16.51) and (16. 9p . Namely, 
has an extra negative sign. 

It is easy to see that the factorized form (16.260 of the Lax operators is 
preserved by the flows of the 2D Toda hierarchy. The Lax equations (16. ip 
can be reduced to the "twisted" version 

{iBC-r),,B-B{{C-'Br).^ 



and 



dtk ' vv ^ />o' 

Qf^ = {{BC-'f)>,C-C{{C-'Bf)^_^, 
(-1)^ = {{CB-'f)^^B- B {{B-'Cf)^^ , 
'-^fgf^ = ((CB-')')<„C - C ((B-'C)')^„ 



^=((^BC-^f\ C-c{{C-^B''^ 



(6.28) 



V / >0 V / >0 

.dB 



lf'd£= ((^CB-^f\ B-B({B-^Cf, 
dtk V / <o V / <o 

.lf'^=((CB-^f\ C-CiiB-^Cf 
dtk V /<o V 



(6.29) 
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of (16 .yp and fl6.12p . These equations, in turn, can be converted to the evolu- 
tion equations 



f-A, (-i)'f -A, (6.30) 

that amount to (16.80 . Thus fl6.26p turns out to define a consistent reduction 
of the 2D Toda hierarchy. The reduced system is substantially the same as 
the Ablowitz-Ladik hierarchy except that the second set of time variables i 
are replaced by L(t) = (— ti, ^2, — ^3, ■ ■ ■)• 

Thus we can conclude that the factorized Lax operators (16.261) with the 
initial data (I6.27P at t = t = persist to be factorized after time evolutions. 
In summary, we have proven the following proposition: 

Theorem 4. The factorization problem Ii6.15\) for the matrix Ii6.20\) yields a 
solution of the 2D Toda hierarchy with the tau function ^-lUj ). The associated 
Lax operators L and L are factorized in the form of 116. 26\} . In particular, 
this is a solution of the Ablowitz-Ladik hierarchy with the second set of time 
variables i being replaced by L{t). 

Remark 3. Unlike the reduction condition (16.51) of Brini et al. [T71 123], our 
reduction condition (16.261) has extra negative signs. This apparent discrep- 
ancy seems to be related to the fact that the usual expression [HI |8] of the 
string amplitude Zx^ of the resolved conifold contains the parameter Q with 
with a negative sign: 

Zx. = J]sA(g-'')5u(g-0(-Q)'"- (6.31) 

Remark 4. The factorization problem (I6.15P for the matrix (I6.19p . too, can 
be treated in the same way. The solution at t = t = reads 

The initial values of the associated Lax operators at t = i = have the 
factorized form 

L = L-^ = e^=(l - q'-^e-^^){l - Qq'-^e-^^). (6.33) 
Such a factorized form is preserved by the fiows of the 2D Toda hierarchy. 
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